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Many radiation shielding problems involve calculations of the response of an isotropic 
detector to radiation of arbitrary angular distribution from uniform rectangular sources. 
In calculations of this type the family of integrals y*s(cos ddS/r 2 )Pi(cosd) and the integral 
J* s (dS /r 2 )exp( — fxt /cosd) are frequently encountered, where 6 is obliquity with respect to an 
axis perpendicular to the plane containing the rectangular radiant surface, S, r is the distance 
from an element of source area, dS, to the detector, n is the attenuation coefficient, and t is 
the barrier thickness. Solutions of the first type of integral facilitate use of Legendre expan- 
sion representations of radiation directional distributions, and may also have application 
in other radiant surface studies, such as illumination and heat exchange engineering. The 
second integral relates to exponentially attenuated radiation from a plane isotropic rec- 
tangular source separated from the detector by a layer of material of thickness t. Formulas, 
expansions, and numerical results are presented. 

1. Introduction 

The penetration of radiation into structures has been the subject of theoretical [1,2] 2 and 
experimental [3] investigations as part of a general program at the National Bureau of Standards 
and other agencies. The theoretical investigations are usually performed in two stages: 
(a) penetration and (b) geometry. The penetration stage refers to generation of data relating 
to interaction of the radiation with barrier materials [4], and geometry refers here to the applica- 
tion of such data to specific configurations of source, barrier, and detector. This paper treats 
the geometrical stage of problems involving rectangular sources, e.g., ceilings of buildings or 
shelters. We draw to some extent from similar treatments in illumination [5 to 8] and heat- 
exchange [9] engineering. The methods developed in this paper, in turn, have application to 
rectangular sources of light and heat as well as gamma radiation. 

In the geometrical stage of many shielding problems, the following situation, generalized 
in figure 1, arises : one lias a radiation detector and a surface S through which radiation emerges. 
In the schematization of reference [1], the flow of radiation through each element of surface, 
dS, is independent of position on the surface, and no scattering of radiation occurs after emer- 
gence from surface S. The intensity and angular distribution, #(cos 0), of radiation crossing 
the plane S is assumed known from the penetration stage of the analysis. The quantity 
</(cos#), usually called the flux [4], is the radiation intensity per solid angle (e.g., photons per 
cm 2 per second per steradian) traveling in the direction 6 with respect to the 2-axis (see fig. 1). 
If the detector has an isotropic response it will register a value I s proportional to an integral 
of #(cos 6) over the total solid angle, Q s , subtended by the surface S from the detector: 

Is=Mdto s (6)g(cose). (1) 

where 

dti s (6)=dS- cos 6/r 2 

and r is the distance from surface element dS to the detector. 

Introducing the quantity ^(0), which is the total range of azimuthal angle, <p, defined by 
the intersection of a 0-cone with the surface S, as in figure 1, we may write (1) in polar coordi- 
nates as 



i Work supported jointly by the Office of Civil and Defense Mobilization, the Navy Bureau of Yards and Docks, and the Defense Atomic 
Support Agency. 

2 Figures in brackets indicate the literature references at the end of this paper. 
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J 1 r+ s (e) 

d(cos0)</(cos0) dip 

= J d(cos6)g(cosd)\l/ s (6). 



(2a) 
(2b) 



As in the example in figure 1, ^(0) may consist of more than one "piece" of azimuth, A i <p s (8), 
so that in the general case it is a summation over all such pieces : 



4> s (8)=A l < Ps (d) + A 2 < Ps (d) + ■ 



(3) 



This quantity is identical with the "angular response function/' ^(0), of reference [1] except for a 
factor 27r. 3 

The surface S may now be specified as a rectangle, as shown in figure 2, with corners at x 
and y coordinates (0, 0), (w, 0), (w, I), and (0, 1). The detector is on the 2-axis, hereafter referred 
to as the "corner position," at z=h. Since the solid angle integration indicated in eq (1) is 
dependent on relative rather than absolute distances we may reduce the number of parameters 
by dividing all dimensions and coordinates by h, i.e., 

a=w/h, b=l/h, a=x/h, and f}=y/h. 

This representation of the relative geometry, rather than representation by horizontal and 
vertical "eccentricities" as in references [1] and [2], is used here in order that the interchange- 
ability of length and width might be used as a partial check of the results. 
Ill Figure 3 shows the use of these scaled coordinates in describing a typical rectangular 
source problem. Radioactive material is uniformly distributed over the roof slab of an under- 
ground shelter, and the problem is to calculate the detector response to radiation transmitted 
by the rectangular ceiling, or rectangular portions of the ceiling, into the enclosure below. 



3 Berger and Lamkin [1] included in \p(fi) a normalizing factor l/2*r, here omitted, 
to their azimuthal range parameter At. 



Thus the quantity \ps(d) defined in (3;, above, is equivalent 
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gure 1. Geometrical schematization of a de- 
tector at distance h from a plane source S of 
arbitrary shape. 



-RADIANT SURFACE, S 




Figure 2. Schematization when S is a rectangle 
of width w and length I with the detector in a 
u corner position" at distance h from the plane 
of S. 
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Figure 3. Underground shelter situation in 
which the coordinates and dimensions in figure 
2 are divided by h, i.e., a = x/h, (3 = y/h, a = w/h 
and b = l/h. 



There are several ways to evaluate the integral (1) for a rectangle, either in polar form 
(2a,b) or using the cartesian parameters a, fi. Berger and Lamkin [1] performed the inte- 
grations numerically, using the polar form, for a few representative angular distributions 
#(cos 0) and a variety of rectangular shapes. Spencer [2] has approximated the function 
^s(0), which may now be written for a rectangle as \l/(d,a,b), by step functions so as to reduce 

(2b) to a combination of values of the cumulative integral d(eos 0) (/(cos 0). 

In this paper, advantage is taken of the fact that r/(cos 0) may usually be represented as a 
sum of Legendre polynomials, i.e., 



g (cos 0)=1b (/+£) QiPi (cos 0), 



where 



(ji= I d(cos 0) g (cos 0) l\ (cos 0) . 
Thus the integral (2b) may be expressed as a sum of elementary integrals 

Ha,b)=±(l + l 2 )(hPi(a,b), 
where the pi(a,b)^ are coefficients of the expansion 



(4) 
(5) 

(6) 
(7) 



The pi(a,b)'s, which are purely geometrical and independent of the penetration stage of the 
problem (i.e., source energy, barrier thickness and material, etc.), may be evaluated as 



Pl (a,b)=\ dQ s (6)P l (cos6) 



(8) 



in which a and b completely specify the surface S. Evaluation of the pt(a,b)'s will be outlined 
in section 4.1. 

The sum (6) converges rapidly if either the Qi or the pi converge rapidly. Since the g/s 
are provided more readily from transport theory [4] than #(cos 0), it is desirable to use (6) 
wherever applicable. 

In the case of penetrating radiations, as from fallout, convergence of (6) maybe inadequate 
as a result of the presence of the unscattered radiation component, which has an angular distri- 
bution of the form 



g°(r y cos 0)=t- exp (— r/cos 0)/cos 0. 

47T 



(9) 
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The superscript ° refers to unscattered radiation, r is the slab thickness, /it, in mean-free-paths, 
and o- is the surface density of radioactive material. Since the g°(T)'& evaluated from (5) and 
(9) converge too poorly for use in the expansion (6) when 0<t~51 } an alternative expansion 

I°(T,a,b)=£-i2T»q n (a,b) (10) 

is useful. The geometry-dependent coefficients q n (a,b), which are obtained by substituting 
(9) into (1), expanding in j n and integrating term by term, are given in section 5.2. 

2. Notation and Definitions 

Notation relevant to this paper, much of which has already been introduced, is listed for 
convenience as follows : 

detector = device with response proportional to the radiation flux at a given point in 
space, 
S= symbolic reference to a specific surface lying in the x,y plane, 
h = distance from surface to detector in the ^-direction, 
w^l— width and length of a rectangular surface, 

x,y= rectangular coordinates parallel to w and I, respectively, with origin at one 
corner of the rectangle, 
corner position = detector position on a line perpendicular to the rectangle at one of its corners, 
i.e., at (0,0,A), 
a=w/h, b=l/h, 
a=x/h, P=y/h, 
dS= element of area of the surface S, 
6= obliquity with respect to the z axis, 

<p= azimuth of the plane through the z axis and the center of dS, measured from 
the x,z plane, 
#s(0)= range of <p covered by all surface elements dS from which the detector can be 
seen at obliquity 0, 
\p(d,a,b) = \f/ s (0) when S refers to a rectangle and the detector is in a corner position, 
r= distance from dS to the detector, 
dQs(0)= solid angle subtended from the detector by element dS at obliquity 
=dS - cos d/r 2 , 
g(cos 0)= radiation flux (e.g., number of photons) traveling in direction per unit solid 
angle, unit time, and unit area normal to direction 0, 
gi= coefficients of a Legendre polynomial expansion of g(cos 6), 
Pi= coefficients of a Legendre polynomial expansion of #s(0)> 
7 lS =response of an isotropic detector to surface S radiating according to #(cos0), 
above, 
I(a J b) = I s when S is a rectangle and the detector is in a corner position, 
I(a,b;a,P) = I s when the detector lies at (ha,hp,h). 
The remaining quantities refer specifically to applications involving penetrating radiations, 
typified by the "fallout problem" shown in figure 3: 

r=iJLt— mean-free-path thickness of layers of material between primary source and the 
electron, where /x is the attenuation coefficient in reciprocal units of the linear 
thickness t, 
<r/4:7r= surface density of contaminant, in the case of the fallout problem, in units of 
photons emitted per unit time from unit surface area into unit solid angle, 
I(T,a,b) = P(T,a,b)+I s (ra,b) 

= the detector response to the unscattered and scattered components, respectively, 
of radiation transmitted by a slab barrier having thickness r and a rectangular 
exit surface of dimensions a, b. 
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3. Arbitrary Detector Positions 



Although the results in this paper relate directly to the response of a detector in a corner 
position, the response at an arbitrary position (a,fi) (shown in fig. 4) is of more general interest. 

A compounding rule found in standard illuminating engineering texts [10] yields the 
detector response, I(a,b;a,/3), at such a position by combining corner-position data from four 
related rectangles : 



I(a,b;a,P) = I(a,p) + I(a-a, b-fi)+I(a, b-p)+I(a-a, 0), 
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Figure 4. Examples using the scaled coordinates of figure 3, of noncorner positions (a) within the rectangular area 

directly under the source and (b) outside this area. 

in which the terms on the right are obtained by replacing the parameters (a,b) of (6) or (10) by 
the quantities indicated. Since both (6) and (10) are odd functions in ab, i.e., 



I(±a,Tb) = -I(a,b), 



(12) 



the rule in (11) applies not only to detector positions 0<o:<a, 0<(3<b as in figure 4a, but to 
any position — oo <«< + oo , — co <^<^-|_ oo as shown in the example in figure 4b. 

4. Legendre Expansion of the Response Function 

4.1. Isotropic Detector 

The primary task performed in this work is the evaluation of the Legendre expansion 
coefficients pi(a,b) by performing the integrations indicated in (8). 
Use of scaled cartesian coordinates leads to the transformations 



and 



cos d=h/r=l/-y/a 2 +l3 2 +l 

dQ s (6)=dS'COse/r 2 =dadf3'(a 2 +(3 2 +l)-* /2 
so that we may write (8) explicitly as 

Pl (a,b)= ( a da fV(a»+j9 I +l)"*.P I (l/Va , +i8 B +l). 

Jo Jo 



(13) 
(14) 

(15) 
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Evaluation of this family of double integrals is outlined in appendix A. The first four pi(a,b)'s 
were found to be 

ab 

p (a,b) = tfm 1 i 2+&2+1 ' ( 16 ) 

i r b .a a , b ~\ 

r>i(ab)=-J tan" 1 , -\ — tan" 1 , y fl7) 

PiW) 2\V6 2 +1 V6 2 +l Va 2 +1 -va 2 +lJ l ^ 

, .. 1 ab f 1 . 1 \ 



and 



Pz(a,b)—-< •—( r tan 1 -=-| — . tan 1 ) 



I ° r ^ 

-Ttz ? tan x . H ^ tan 1 _ 

2 L(& 2 +i) f V& 2 +i (a 2 +i)^ V^+i 



5Tb , a . a _, 6 

I _ ran x 1 ton - 1 

hi) 
ab 



a 2 +6 2 +l 



C-^I+FTi)]}- (19) 



The po(a,b) term is the result given by Spencer [2] for a diffuse (i.e., g(cos 0) = constant) rec- 
tangular source and isotropic detector, and pi(a,b) is the result given by illumination [5-8] and 
heat-exchange [9] engineering studies for a diffuse rectangular source and a surface-type detector 
parallel to the source plane. 

The terms become progressively more complicated with increasing I, but the even and odd 
general terms as given by (38) and (42) in appendix A are amenable to programing for an 
automatic computer. Table 1 contains pi(a,bys evaluated using an IBM type 704 computer 
for 0<Z<13 over a grid of a and b values 0.1 <a<6<20. These may now be inserted in (6) to 
give the detector response I(a,b) in a corner position. 

These pz(a,6)'s may be checked by reconstructing \[/(6,a,b) according to (7) and comparing 
with the exact results given in reference [1]. As shown in figure 5, ^ (0,a,b) is reasonably well 
described for the case a=l, b = 2, neglecting terms in excess of 1=7. If, however, higher g x 
terms are available and contain appreciable information, as for a severely peaked angular 
distribution, greater accuracy may be achieved by including additional terms in (6). 

4.2. Cosine Detector 

The cji coefficients may in principle be generated for any azimuthally independent uniform 
plane source by performing the integration indicated in (5) by analytic, numerical, or graphical 
means. Thus the Legendre expansion method might also be fruitful in similar rectangular 
source studies in illumination and heat-exchange engineering in which the g(cos 6) is known. 
In these studies, however, the radiation field is usually described in terms of a surface-element 
detector ("plane detector") whose response is proportional to the cosine of the angle of in- 
cidence. 

The orientation of such a plane detector may be specified by its direction cosines A, B, 
and C with respect to the a, fi, and z axes. Legendre expansion coefficients pi,A,B,c(a,b) 
(analogous to the isotropic detector pi(a } b)'s) may each be written as the sum of three com- 
ponents (see ref [8], eq (17)), 

Pi,A,B,c(a,b)=Ap la (a,b) + Bp l( 3(a,b) + Cp l2 (a,b), (20) 

where p la (a,b) and p# (a,b) are the Legendre coefficients for a plane detector perpendicular 
to S and to the a- or 0-axis, respectively, and the p iz (a,b)'$ are the coefficients for a detector 
parallel to S, i. e., perpendicular to the 2-axis. 
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Table 1.* Coefficients p/(a, b) of the Legendre expansion (7) of the angular response function \p(d, a, b) 



I 


6 = 0.1 


6=0.2 


6=0.2 


6=0.5 


6=0.5 


6=0.5 




a=0.1 


0=0.1 


a=0.2 


a=0.1 


a =0.2 


a=0.5 





9.9012 (-3) 


1.9516 (-2) 


3.8471 (-2) 


4.4514 (-2) 


8.7819 (-2) 


2.0136 (-1) 


1 


9. 8685 


1.9358 


3. 7980 


4.2915 


8. 4281 


1 . 8807 


2 


9. 8034 


1.9046 


3. 7010 


3. 9869 


7. 7547 


1 . 6330 


3 


9. 7064 


1 . S5S5 


3. 5587 


3. 5656) 


6. 8256 


1.3035 


4 


9. 5781 


1.7984 


3. 3750 


3. 0658 


5. 7261 


9.3443 (-2) 


5 


9. 4193 


1 . 7254 


3. 1545 


2. 5301 


4. 5524 


5. 7010 


6 


9. 2312 


L.6408 


2. 9031 


2.0015 


3. 3997 


2. 4985 


7 


9. 01 51 


1.5461 


2. 6270 


1.5177 


2. 3514 


2.2553 (-4) 


8 


8. 7724 


1.4430 


2. 3332 


1.1076 


1. 4701 


-1.5825 (-2) 


9 


8. 5047 


1 . 3333 


2. 0287 


7.8884 (-3) 


7.9220 (-3) 


-2. 3202 


10 


8.2140 


1.2188 


1.7206 


5. 6642 


3. 2668 


-2.3214 


11 


7. 9020 


1.1015 


1.4159 


4. 3403 


5.7140 (-4) 


-1.8030 


12 


7.5710 


9.-6316 (-3) 


1.1211 


3. 7609 


-5.2395 


-1.0173 


13 


7. 2231 


8. 6565 


8.4197 (-3) 


3. 7105 


-4.9818 


-2.0192 (-3) 



I 


6 = 1.0 


6 = 1.0 


6 = 1.0 


6 = 1.0 


6=2.0 


6=2.0 




a=0.1 


a=0.2 


a=0.5 


a-1.0 


a=0.1 


a=0.2 





7.0418 (-2) 


1.3912 (-1) 


3.2175 (-1) 


5.2360 (-1) 


8.9117 (-2) 


1.7632 (-1) 


1 


6.3910 


1.2572 


2. 8332 


4. 3521 


7. 4970 


1.4769 


2 


5. 2552 


1.0233 


2. 1667 


2. 8868 


5. 3170 


1.0348 


3 


3. 9093 


7.4574 (-2) 


1.3858 


1.3137 


3. 2627 


6.1652 (-2) 


4 


2. 6458 


4. 8456 


6.6728 (-2) 


8.0187 (-3) 


1.9494 


3.4619 


5 


1.6848 


2. 8473 


1.4080 


-5.6244 (-2) 


1 . 4458 


2.3811 


6 


1.1172 


1.6480 


-1.4393 


-6. 2546 


1.3947 


2. 2105 


7 


8.9990 (-3) 


1. 16(11 


-2. 1736 


-3. 2689 


1.3716 


2.1019 


8 


8. 9914 


1.1114 


-1.6157 


4.4199 (-3) 


1.1807 


1 . 6654 


9 


9. 5780 


1. IS50 


-6.8461 (-3) 


2.6786 (-2) 


8.9548 (-3) 


1.0498 


10 


9. 5857 


1.1468 


-3.5674 (-4) 


2.7119 


6. 8459 


5.9365 (-3) 


11 


8. 5803 


!). 1332 (-3) 


1.0787 (-3) 


1.1910 


6.2105 


4. 4251 


12 


6. 8321 


5. 4331 


-1.0131 


-5.9052 (-3) 


6. 3528 


4.5618 


13 


4.9912 


1.6968 


-3. 4306 


-1.5587 (-2) 


6. 1304 


4. 0585 



/ 


6=2.0 


6 = 2.0 


6=2.0 


6 = 5.0 


6 = 5.0 


6=5.0 




a=0.5 


o=1.0 


a=2.0 


a=0.1 


a=0.2 


a=0.5 





4.1152 (-1) 


6.8472 (-1) 


9.2730 (-1) 


9.7727 (-2) 


1.9351 (-1) 


4.5396 (-1) 


1 


3. 3564 


5. 2582 


6. 5269 


7. 7896 


1. 5352 


3. 4997 


2 


2. 1822 


2. 8577 


2. 6667 


5. 0420 


9.7983 (-2) 


2.0456 


3 


1 . 0654 


7.0518 (-2) 


-2.6030 (-2) 


2.9189 


5. 4792 


8.9638 (-2) 


4 


3.3703 (-2) 


-4. 0123 


-1.2049 (-1) 


1.9463 


3. 4566 


3.3712 


5 


4.3254 (-3) 


-0. 6172 


-6.0731 (-2) 


1.7051 


2. 8992 


1.7192 


6 


7.2294 (-4) 


—2. 9200 


2.8998 


1.5548 


2. 5295 


8.4961 (-3) 


7 


1.3615 (-3) 


8.1359 (-3) 


5. 6423 


1.2639 


l.ss.V.t 


-4. 1356 


8 


-3. 9207 


1.9132 (-2) 


1.9243 


9.8612 (-3) 


1.2767 


-1.3530 (-2) 


9 


-1.1888 (-2) 


9.8109 (-3) 


-2. 5299 


8.6130 


9.8248 (-3) 


-1. 3389 


10 


-1.4791 


-8.7945 (-4) 


-3. 1339 


8. 1909 


8. 6299 


-8.0041 (-3) 


11 


-9.9201 (-3) 


-2.6009 (-3) 


-3.7923 (-3) 


7. 3200 


6. 6370 


-4. 5171 


12 


-1.6638 


9.0735 (-4) 


2.0619 (-2) 


5. 9839 


3.8145 


-3. 6896 


13 


3. 7038 


1.4956 (-3) 


1.7773 


5. 0339 


1.8646 


-1. 7941 



I 


6 = 5.0 


6 = 5.0 


6=5.0 


6 = 10.0 


6 = 10.0 


6=10.0 




a=1.0 


a=2.0 


a = 5.0 


a=0.1 


a=0.2 


a=0.5 





7.6616 (-1) 


1.0697 (0) 


1.2925 (0) 


9.9172 (-2) 


1.9640 (-1) 


4.6117 (-1) 


1 


5. 5285 


6.9768 (-1) 


7.6053 (-1) 


7.8117 


1.5396 


3.5108 


2 


2. 5907 


2. 1768 


1.3464 


4. 9749 


9.6)643 (-2) 


L'. (11 21' 


3 


3.8325 (-2) 


-8.0873 (-2) 


-1. 60C7 


2. 8872 


5.4156 


8.8055 (-2) 


4 


-4. 7545 


-1.1485 (-1) 


-9.3400 (-2) 


1.9878 


3. 5396 


3. 5782 


5 


-4. 1089 


-1.5559 (-2) 


5. 5039 


1.7418 


2. 9725 


1. 9018 


6 


-1.4967 


5. 0088 


6. 8498 


1.5308 


2.4815 


7.2960 (-3) 


7 


-3.3765 (-3) 


3.2101 


-1.4585 


1.2260 


1.8101 


- 6. 0253 


8 


6.9961 (-4) 


-1. 1825 


-4. 9763 


9.9606 (-3) 


1.2966 


-1.3032 (-2) 


9 


7.9552 (-3) 


-2. 3266 


-4.7962 (-3) 


8. 9743 


1.0547 


-1.1588 


10 


1.2904 (-2) 


-4.7751 (-3) 


3.4526 (-2) 


8. 2073 


8.6626 (-3) 


-7.9255 (-3) 


11 


7.2894 (-3) 


9. 7856 


1.4253 


7. 0017 


6. 0006 


-6. 1052 


12 


-4. 0976 


6. 8161 


-2. 1923 


5. 8764 


3. 5995 


-4. 2231 


13 


-9. 4124 


-1.0631 


-1.8035 


5. 2905 


2. 3776 


-5.1238 (-4) 


I 


6 = 10.0 


6 = 10.0 


6 = 10.0 


6 = 10.0 


6=20.0 


6 = 20.0 




0=1.0 


a=2.0 


a = 5.0 


a = 10.0 


a=0.1 


a=0.2 





7.8045 (-1) 


1.0973 (0) 


1.3500 (0) 


1.4300 (0) 


9.9544 (-2) 


1.9715 (-1) 


1 


5. 5503 


7.0183 (-1) 


7.6864 (-1) 


7.7903 (-1) 


7.8146 


1. 5402 


2 


2. 5244 


2. 0484 


1.0771 


6.9836 (-2) 


4. 9567 


9.6278 (-2) 


3 


3.5197 (-2) 


-8.6847 (-2) 


-1.7179 


-1.8686 (-1) 


2. 8829 


5. 4071 


4 


-4. 3432 


-1.0684 (-1) 


-7.6187 (-2) 


-5.1368 (-2) 


2. 0009 


3. 5657 


5 


-3. 7478 


-8.6468 (-3) 


6. 8721 


8. 6458 


1. 7470 


2. 9830 


6 


-1. 7360 


4.5362 (-2) 


5. 7789 


4.1505 


1. 5206 


2. 4612 


7 


-7.1162 (-3) 


2. 4916 


-2. 9C08 


-4. 7938 


1.2200 


1. 7982 


8 


1. 7058 


-9.7509 (-3) 


-4. 4287 


-3. 4798 


1.C040 


1.3126 


9 


1.1525 (-2) 


-1.6370 (-2) 


9.3314 (-3) 


2. 8212 


9.0385 (-3) 


1.0675 


10 


1. 3040 


-4.6562 (-3) 


3.3426 (-2) 


2. 9638 


8. 1452 


8.5385 (-3) 


11 


4.1349 (-3) 


3. 6459 


1.2963 (-3) 


-1.6478 


6. 9345 


5. 8661 


12 


-5. 1356 


. 4.9501 


-2.4409 (-2) 


-2. 5371 


5. 9230 


3. 6927 


13 


-6. 8582 


3.9658 


-6.9441 (-3) 


8.8627 (-3) 


5. 3592 


2.5151 



See footnote at end of table, p. 128. 
5398S9— 60 4 
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Table 1.* Coefficients Pi(sl, b) of the Legendre expansion (7) of the angular response function 

rp(d, a, b) — Continued 



I 


6=20.0 


6=20.0 


6=20.0 


6=20.0 


6=20.0 


6=20.0 




a=0.5 


a=1.0 


a=2.0 


a = 5.0 


a = 10.0 


a =20.0 





4.6302 (-1) 


7.8415 (-1) 


1.1047 (0) 


1.3673 (0) 


1.4593 (0) 


1.5002 (0) 


1 


3. 5122 


5. 5532 


7.0240 (-1) 


7.6994 (-1) 


7.8112 (-1) 


7.8379 (-1) 


2 


2. 0031 


2. 5062 


2. 0124 


9.9215 (-2) 


5.5376 (-2) 


3.5245 (-2) 


3 


8.7841 (-2) 


3.4772 (-2) 


-8.7686 (-2) 


-1.7372 (-1) 


-1.8997 (-1) 


-1.9395 (-1) 


4 


3. 6436 


-4. 2128 


-1.0426 (-1) 


-7.0077 (-2) 


-4.0924 (-2) 


-2.6306 (-2) 


5 


1. 9280 


-3. 6957 


-7.6179 (-3) 


7. 1094 


9. 0285 


9. 5183 


6 


6.7898 (-3) 


-1. 8370 


4.3358 (-2) 


5. 3034 


3.3318 


2. 1754 


7 


-6. 3214 


-7.7066 (-3) 


2. 3750 


-3. 1702 


-5.2295 


-5. 7885 


8 


-1.2633 (-2) 


2. 5016 


-8.1707 (-3) 


-4.0518 


-2. 8235 


-1. 8836 


9 


-1. 1268 


1.2164 (-2) 


-1.5107 (-2) 


1.2254 


3. 2956 


3. 9067 


10 


-8.2358 (-3) 


1. 2420 


-5.8884 (-3) 


3. 0463 


2. 4392 


1.6728 


11 


-6. 4412 


3.4648 (-3) 


2. 3212 


-1.7760 (-3) 


-2. 1487 


-2. 7972 


12 


-3. 9897 


-4. 6694 


5. 8792 


-2.2148 (-2) 


-2.1264 


-1.5087 


13 


-1.6924 (-4) 


-6. 1736 


5. 3201 


-3.7939 (-3) 


1.4029 


2. 0756 



*Note: The figures in parenthesis in tables 1,2, 3, and 4 indicate the power of ten by which the numbers alongside and below in the same 
column are to be multiplied; e.g., 9.9012(-3) =0.0099012. 




Figure 5. Reconstruction of the response 
f unction \p(d, a, b) according to (10) for a=b=l 
using terms 0< /< 7 (dotted line) and 0</<9 
(dashed line). The exact function as given by 
Berger and Lamkin [1] is shown as the solid 
line. 



cos 9 



The dependence on the angle of incidence of radiation on a plane detector for geometries 
described by p la (a,b) and p# (a,b) reduces to a cos 9 and cos 0, respectively, modifying the 
integrand in (15) to give 



and 



Pia(a, &) = f dfi ^ ^aP,(l/Va 2 +)S 2 +l)/(a 2 +^+l) 2 
Jo Jo 

P*(a, b) = f " da P pd(3Pi(lN* 2 +P 2 +l)Ka 2 +F+iy. 
Jo Jo 



(21) 



(22) 



Evaluation of these integrals should be similar to that for (15). In the geometry described 
by the p Xz (a,6)'s, the angles of emission and incidence are identical, so that the integrand in 
(15) requires only an additional factor of cos 0. Thus, using the relation 



cos dP t (cos 0) = {«>,_! (cos 0) + (Z+l)P z+1 (cos 0)}/(2Z+l), 
the p u (a,b)'s may be obtained from the p z (a,b)'s in table 1 as 

P l2 (a,b) = {lp l . 1 (a,b) + (l + l)p l+1 (a } b)}/(2l + l). 
5. Unscattered Radiation Component 



(23) 



(24) 



In cases involving penetrating radiations, such as gamma rays from the plane source 
shown in figure 3, it is usually desirable to make a separate calculation of the unscattered 
radiation component, P(T,a,b). Two reasons for this separation of components are: 
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(a) In contrast to the scattered radiation component, P(r,a,b), the unscattered com- 
ponent, I°(T,a,b), is analytic in nature and independent of the barrier material, so that it may- 
be calculated once and for all for the relevant parameters r, a, and b. 

(b) Since the angular distribution, g°(r, cos 0), of unscattered radiation from a plane 

isotropic source as expressed in (9) tends toward j- (1/cos 6) as r becomes small, representation 



as a Legendre expansion in positive powers of cos 6 becomes impractical, 
usually does not arise for the scattered component. 

5.1. t>1: Legendre Expansion 



This difficulty 



In order to examine the behavior discussed in (b), above, #z°( r )' s were evaluated by sub- 
stituting (9) into (5). The result 



g l °{r)=^{P l E ]+l {r)} 



where Ej(t) is the exponential integral [11] 



Ej(r)- 



-s; 



e~ Tfi u~ j du, 



(25) 



(20) 



and P^Ej+^t)} is obtained by replacing (cos d) j inP t (cos 6) by E j+1 (r). Since these g t ° (t), 
although of fairly general interest, do not seem to appear explicitly in the literature, a tabulation 

of the quantity— (/ + |)^°(t), (i.e., (Z+J) P t {E j+l (r)}), for <l< 13 and 0.001 < r < 100 is 
a* 

presented in table 2. 

A check on the values in table 2 is shown in figure 6, in which (47r/o-) e T g° (r, cos 6), recon- 
structed according to (4) for 0<l<7 and 0.2 <r< 10, are indicated by the dotted curves. 
Comparison may now be made with similar results from the exact expression (9), shown by the 

Table 2. Coefficients (47r/<r)(Z + i)gJ(r) of a Legendre expansion (4) of the angular distribution (47r/o-)g' > (r, cos 0) 

of unscattered radiation from a plane isotropic source 

The factor iir/a implies unit source strength per steradian per unit of source area, and (l-\- l /2) is included to allow cumulative multiplication 
directly with the pi {a,b)'s (see table 1) according to (6). 



\r 


0.001 


0.01 


0.1 


0.2 


0.5 


1.0 


_1\ 

















3.1943 (0) 


2.0190 (0) 


9.1145 (-1) 


6.1135 (-1) 


2.7989 (-1) 


1.0970 (-1) 


1 


1. 4889 


1. 4245 


1.0838 (0) 


8. 6130 


4. 8996 


2. 2275 


2 


-0.1143 


-3. 2090 


-7.1755 (-1) 


-2. 0852 


1. 3130 


1. 3712 


3 


-2. 2989 


-2. 1124 


-1.2757 (0) 


-8. 3188 


-2. 6901 


-2.6555 (-2) 


4 


7. 2752 


3. 3973 


3.6314 (-1) 


-9.1715 (-2) 


-2. 1634 


-9. 3776 


5 


2. 8703 


2. 5450 


1.2524 (0) 


6.7606 (-1) 


1.0298 


-3. 3952 


6 


-8. 0184 


-3. 4034 


-9.6551 (-2) 


2. 4954 


1. 9507 


4. 0965 


7 


-3. 3335 


-2. 8604 


-1.1624 (0) 


-5. 1631 


-8.9640 (-3) 


3. 9813 


8 


8. 5595 


3. 3377 


-9.7963 (-2) 


-3. 2542 


-1.4873 (-1) 


-9.6373 (-3) 


9 


3. 7303 


3. 1040 


1.0492 (0) 


3. 7667 


-3.8228 (-2) 


-3.0111 (-2) 


10 


-8. 9790 


-3. 2372 


2.3859 (-1) 


3. 5348 


1.0293 (-1) 


-4.9996 (-3) 


11 


-4. 0810 


-3. 2983 


-9.3042 


-2. 6113 


5.7819 (-2) 


1.8731 (-2) 


12 


9. 3161 


3. 1186 


-3.3868 


-3. 5306 


-6. 4944 


1.0139 


13 


4. 3971 


3. 4559 


8. 1428 


3. 3172 


-6. 1991 


-9.7479 (-3) 


\ 

\ T 


2.0 


5.0 


10.0 


15.0 


20.0 


100.0 


'\ 

















2.4451 (-2) 


5.7415 (-4) 


2.0785 (-6) 


9.5930 (-9) 


4.9178 (-11) 


1.8418 (-46) 


1 


5. 6301 


1.4947 (-3) 


5. 7483 


2.7162 (-8) 


1.4107 (-10) 


5. 4717 


2 


5. 1875 


1. 8564 


8.1118 


4. 0293 


2. 1490 


8. 9433 


3 


2. 1895 


1. 6196 


8. 8022 


4. 7268 


2. 6263 


1.2159 (-45) 


4 


-6.1965 (-3) 


1.0195 


7. 9587 


4. 7804 


2.8115 


1. 5034 


5 


-1.4221 (-2) 


3.9434 (-4) 


6. 1188 


4. 2970 


2. 7243 


1. 7498 


6 


-5.7195 (-3) 


-1.1402 (-5) 


3. 9668 


3. 4702 


2. 4242 


1. 9500 


7 


3. 9476 


-1.4339 (-4) 


2.0817 


2. 5187 


1. 9934 


2. 1009 


8 


5. If06 


-9.8609 (-5) 


9.1367 (-7) 


1. 6300 


1. 5167 


2. 2013 


9 


5.1102 (-4) 


-1.3120 


8.0998 (-8) 


9.1882 (-9) 


1. 0653 


2. 2522 


10 


-2. 6394 


3. 0530 


-1.5133 (-7) 


4. 3555 


6.8479 (-11) 


2.2561 


11 


-1. 5451 


2. 4977 


-1.4134 


1. 4721 


4. 0174 


2.2172 


12 


8. 7741 


3.8280 (-6) 


-5.0961 (-8) 


1.6851 (-10) 


2. 0205 


2. 1418 


13 


1.2960 (-3) 


-9.5729 (-6) 


4.5106 (-9) 


-2.5223 (-11) 


9.8335 (-12) 


2. 0335 
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Figure 6. The angular distribution of unscat- 
tered radiation from a plane isotropic source 
reconstructed according to (4) and (25) (dotted 
lines) as compared with the exact expression 
(9) (solid lines) for 0.2<r<10. 



cos 8 

solid curves. From both table 2 and figure 6 it can be seen that the Legendre expansion 
approach, using tables 1 and 2, should give reliable results for unscattered as well as scattered 
radiation in the region 1^7^15. 

5.2. r"<l: Power Series Expansion 

The parameter r may in some cases be very small, and for this reason it seemed desirable 
to find a method for evaluating I°(T,a,b), i.e., 



F(T,a,b)=~ j (dS/f*) exp (-r/cos 6), 



(27) 



for small r. No success was had in obtaining a closed-form result, but q n (a,b) coefficients of 
the r n power series in (10) were evaluated (see appendix B) which provide adequate convergence 
of (10) from r=0 up to values of r for which the g° t (a,b)'s converge rapidly. 
The first five q n (a,b) terms are 

T . b ^ 1 / a 2 Y +1 /V<?+!* iV^+1 fc* 2 2 >(j!) 2 / a 2 +l V +1 \ 

(28) 
in which the inner sum overj vanishes for i = 0, 



b sinh l -\-a sinh" 1 



y& 2 +i 



V^+i 



-tan 



ab 



Va 2 + 6 2 + 



i}- 



&(«,&)= 



ab 
W 



(29) 
(30) 



w)-4{-*#^[w+« sb h-.^ + * + .) ,^^,^imi}, 



2s(ffl,6)= 



(31) 
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and 



24(o > 6)=jj/(2!) 2 ,(o,6)+ 



"3-l" 



1.3/' 



(32) 



For large a and 6, (a<b), the expression within the braces in (28) becomes small very 
quickly with increasing i. Indeed, over the entire range of a and b, i.e., <a<&<°°, summation 
over i from i=0 to f=6 gives values for g {dfi) departing less than 0.01 percent from their 
respective values when the summation is carried to i\>20. A partial check of (28) is obtained 
by noting that the series part of the expression vanishes for large b, that is, 



lim q fl (a,b)=-sm]\ 1 a> 

6->oo * 



(33) 



in agreement with the result of Andrews [12] for the detector response to an infinite strip source 
in a nonattenua ting medium. 

The remaining q n (a,b) terms form two distinct sets, as did the j)i (a,&)'s, depending on 
whether n is even or odd. The general terms for these sets are given in (57) and (58) of 
appendix B. 

Values of 1° (r,a,6)4ir/<r for 0.01 <r< 10 and 0.1 <a<b< 20, as evaluated with the help 
of an IBM 704 computer, are presented in table 4. Numbers preceded by asterisks were 
evaluated using the Legendre expansion approach (15) and (25) for 0</< 13, and the others 
using the r n series (10) for ()<// < 20. 

Since evaluation of (10) for small specific values of r may be faster than interpolation from 
table 4, q n (a,b) is given in table 3 for <n<9 and 0.1 <a<6< 20. 



Table 3. Coefficients q n (a,b) of the r n power series expansion (10) of the corner position flux 1° (r,a,b) of un- 
scattered radiation from a rectangular plane isotropic source 



11 


6=0.1 


6=0.2 


6=0.2 


6 = 0.5 


6 = 0.5 


6 = 0.5 




a=0.1 


fl=0.1 


a=0.2 


a =0.1 


a=0.2 


a = 0.5 





9.9340 (-3) 


1.9675 (-2) 


3.8971 (-2) 


4.0222 (-2) 


9. 1003 (-2) 


2. 1600 (-1) 


1 


-9. 9009 


-1.9830 


-3. 9481 


-4. 8047 


-9. 5050 


-2.3215 


2 


5. ()()()() 


10. 0000 (-3) 


2. 0000 


2. 5000 


5. 0000 


1.2500 


3 


-1. 0722 


-3.3609 


-0.7547 (-3) 


-8.0819 (-3) 


-1.7443 


-4.4900 (-2) 


4 


4. 1944 (-4) 


8.4722 (-4) 


1.7111 


2. 2039 


4.5094 (-3) 


1.2153 


5 


-8.4109 (-5) 


- 1 . 7080 


-3.4679 (-4) 


-4.7282 (-4) 


-9.5873 (-4) 


-2.0333 (-3) 


6 


1.4075 


2.8710 (-5) 


5.8574 (-5) 


8.2390 (-5) 


1.0783 


4.7047 (-4) 


7 


-2.0174 (-0) 


-4.1308 (-6) 


-8.4805 (-0) 


-1. 2320 


-2.5213 (-5) 


-7.4048 (-5) 


8 


2.5302 (-7) 


5.2148 (-7) 


1. 0744 


1.6141 (-0) 


3.3183 (-6) 


1.0090 


9 


-2.8208 (-8) 


-5.8435 (-8) 


-1.2101 (-7) 


-1.8819 (-7) 


-3.8869 (-7) 


-1.2248 (-6) 


n 


6=1.0 


6=1.0 


6=1.0 


6 = 1.0 


6 = 2.0 


6 = 2.0 




a=0.1 


a=0.2 


a=0.5 


a = 1.0 


o = 0.1 


a=0.2 





7.8327 (-2) 


1.5540 (-1) 


3.0887 (-1) 


0.3951 (-1) 


1.1040 (-1) 


2.1940 (-1) 


1 


-8. 8020 


-1.7535 


-4. 2718 


-7. 9336 


-1.4421 


-2.8755 


2 


5. 0000 


10.0000 (-2) 


2. 5000 


5. 0000 


10.0000 (-2) 


2. 0000 


3 


-1.9154 


-3. 8455 


-9.8024 (-2) 


-2. 1346 


-4. 9338 


-9.8916 (-2) 


4 


5.5094 (-3) 


1. 1222 


2. 9514 


6.9444 (-2) 


1. 9472 


3.9111 


5 


-1.3114 


-2.0517 (-3) 


-7.1478 (-3) 


-1.8357 


-0.5195 (-3) 


-1.3113 


6 


2.0050 (-4) 


5.2848 (-4) 


1.4593 


4.1049 (-3) 


1.9117 


3.8496 (-3) 


7 


-4.4892 (-5) 


-9. 1365 (-5) 


-2.5829 (-4) 


-7.9826 (-4) 


-5.0012 (-4) 


-1.0081 


8 


0.8492 (-0) 


1. 3982 


4.0448 (-5) 


1. 3771 


1.1811 


2.3829 (-4) 


9 


-9.3939 (-7) 


-1.9233 (-6) 


-5.0906 (-0) 


-2.1392 (-5) 


-2.5399 (-5) 


-5.1290 (-5) 


ii 


6 = 2.0 


6 = 2.0 


6 = 2.0 


6 = 5.0 


6 = 5.0 


6 = 5.0 




a =0.5 


c = 1.0 


a = 2.0 


a = 0.1 


a=0.2 


a=0.5 





5.2530 (-1) 


9.3072 (-1) 


1.4074 (0) 


1.3708 (-1) 


2.7203 (-1) 


6.5729 (-1) 


1 


-7.0450 


-1.3285 (0) 


-2.2916 


-2.3108 


-4.0120 


-1.1372 (0) 


2 


5. 0000 


1.0000 


2. 0000 


2. 5.000 


5. 0000 


1. 2500 


3 


-2.5139 


-5.3007 (-1) 


-1.2410 


-2.3179 


-4. 6397 


-1. 1666 


4 


1.0069 


2. 2222 


0.1111 (-1) 


1.9451 


3. 8944 


9.8090 (-1) 


5 


-3.4104 (-2) 


-7.8243 (-2) 


-2.5253 


-1.4085 


-2.9404 


-7. 4123 


6 


1. 0095 


2. 3951 


9.0494 (-2) 


9.9117 (-2) 


1. 9849 


5. 0081 


7 


-2.6025 (-3) 


-6.5116 (-3) 


-2.8734 


-5.9935 


-1.2004 


-3. 0320 


8 


0.3360 (-4) 


1. 5949 


8.2096 (-3) 


3. 2071 


6.5448 (-2) 


1. 6550 


9 


-1.3730 


-3.5550 (-4) 


-2. 1351 


-1.0173 


-3. 2405 


-8.2045 (-2) 
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Table 3. Coefficients q n (a,b) of the r n power series expansion (10) of the corner position flux 1° (r,a,b) of un- 
scattered radiation from a rectangular plane isotropic source — Continued 



n 


6 = 5.0 


6 = 5.0 


6=5.0 


6 = 10.0 


6 = 10.0 


6 = 10.0 




a = 1.0 


a = 2.0 


a = 5.0 


a= 0.1 


a= 0.2 


a= 0.5 





1.1879 (0) 


1.8793 (0) 


2.7262 (0) 


1.4685 (-1) 


2.9217 (-1) 


7.0607 (-1) 


1 


-2. 1837 


-3. 9333 


-7. 3832 


-2. 9966 


-5. 9833 


-1.4798 (0) 


2 


2. 5000 


5. 0000 


1.2500 (1) 


5. 0000 


10. 0000 


2. 5000 


3 


-2.3782 


-5.0728 


-1. 6595 


-8. 6256 


-1.7256 (0) 


-4. 3227 


4 


2. 0139 


4. 4444 


1. 8403 


1.4307 (0) 


2. 8622 


7. 1701 


5 


-1.5272 


-3. 4397 


-1. 7554 


-2. 1472 


-4. 2957 


-1.0762 (1) 


6 


1. 0353 


2.3716 


1. 4695 


2. 8721 


5. 7461 


1. 4399 


7 


-6.2923 (-1) 


-1.4665 


-1.0971 


-3. 4332 


-6. 8689 


-1. 7217 


8 


3. 4493 


8.1909 (-1) 


7.3983 (0) 


3. 6949 


7. 3929 


1. 8536 


9 


-1.7178 


-4. 1614 


-4. 5541 


-3. 6090 


-7. 2214 


-1.8112 


n 


6 = 10.0 


6 = 10.0 


6 = 10.0 


6 = 10.0 


6 = 20.0 


6 = 20.0 




a= 1.0 


a= 2.0 


a= 5.0 


a = 10.0 


a= 0.1 


a= 0.2 





1.2849 (0) 


2.0692 (0) 


3. 1467 (0) 


3.7961 (0) 


1.5182 (-1) 


3.0211 (-1) 


1 


-2. 8671 


-5. 2864 


-1.0571 (1) 


-1.6127 (1) 


-3. 6878 


-7.3658 


2 


5. 0000 


1.0000 (1) 


2. 5000 


5. 0000 


10. 0000 


2.0000 (0) 


3 


-8. 7046 


-1.7838 


-5. 0356 


-1.2892 (2) 


-3.3683 (0) 


-6. 7373 


4 


1.4444 (1) 


2. 9722 


8. 8889 


2. 8194 


1.1195 (1) 


2.2391 (1) 


5 


-2. 1687 


-4.4719 


-1.3854 (2) 


-5.3227 


-3. 3586 


-6. 7176 


6 


2. 9038 


6. 0067 


1. 9201 


8. 8285 


8. 9635 


1.7928 (2) 


7 


-3. 4758 


-7. 2190 


-2.3848 


-1.3066 (3) 


-2.1364 (2) 


-4. 2732 


8 


3. 7462 


7. 8159 


2. 6756 


1. 7475 


4. 5831 


9. 1672 


9 


-3. 6647 


-7. 6824 


-2. 7323 


-2. 1341 


-8.9221 


-1.7846 (3) 


n 


6 = 20.0 


6 = 20.0 


6 = 20.0 


6 = 20.0 


6 = 20.0 


6 = 20.0 




a= 0.5 


a= 1.0 


a= 2.0 


a= 5.0 


a = 10.0 


« = 20.0 





7.3091 (-1) 


1.3345 (0) 


2.1679 (0) 


3.3843 (0) 


4.2228 (0) 


4.8801 (0) 


1 


-1.8253 (0) 


-3. 5577 


-6. 6641 


-1.3955 (1) 


-2.2546 (1) 


-3.3719 (1) 


2 


5. 0000 


1.0000 (1) 


2.0000 (1) 


5. 0000 


1.0000 (2) 


2.0000 (2) 


3 


-1.6854 (1) 


-3.3781 


-6. 8107 


-1.7800 (2) 


-3.9741 


-1.0231 (3) 


4 


5. 6007 


1.1222 (2) 


2.2611 (2) 


5. 9444 


1.3972 (3) 


4. 4611 


5 


-1.6803 (2) 


-3. 3669 


-6. 7850 


-1.7901 (3) 


-4. 3433 


-1.6801 (4) 


6 


4. 4847 


8. 9882 


1.8127 (3) 


4. 8099 


1.2006 (4) 


5. 5605 


7 


-1.0690 (3) 


-2.1430 (3) 


-4. 3264 


-1.1564 (4) 


-2.9741 


-1.6422 (5) 


8 


2. 2935 


4. 5990 


9. 2953 


2. 5046 


6. 6542 


4. 3835 


9 


-4. 4652 


-8. 9564 


-1.8123 (4) 


-4. 9253 


-1.3550 (5) 


-1.0685 (6) 



Table 4. Unscattered flux (4:Tr/a)I°(r,a,b) of radiation, at a corner position, from a rectangular plane isotropic 
source calculated using tables 1 and 2 (numbers preceded by asterisks) and using table 3 (no asterisks) 



T 


6=0.1 
a =0.1 


6=0.2 
a=0.1 


6 = 0.2 
a=0.2 


6=0.5 
a=0.1 


6 = 0.5 
a=0.2 


6=0.5 
a=0.5 




0.01 
0. 02 
0.05 
0.1 


9.93 (-3) 

9.83 

9.74 

9.44 

8.99 


1.97 (-2) 

1.95 

1.93 

1.87 
1.78 


3.90 (-2) 

3.86 

3.82 

3.70 

3.52 


4.62 (-2) 

4.57 

4.53 

4.39 

4.17 


9.16 (-2) 

9.07 

8.97 

8.69 

8.25 


2.16 (-1) 

2.14 

2.11 

2.05 

1.94 


0.2 


8.13 
*8.48 


1.61 
*1.65 


3.18 
*3.21 


3.75 
*3.78 


7.43 
*7.44 


1.74 
•1.75 


0.5 


6.02 
*5.63 


1.19 
•1.15 


2.35 
*2.31 


2.75 

*2.72 


5.44 
*5.44 


1.26 
*1.26 


1.0 


3.64 
*3.66 


7.18 (-3) 
*7.02 


1.42 (-2) 
•1.42 


1.64 
*1. 64 


3.23 
*3.23 


7.38 (-2) 
*7.38 


2.0 


1.34 
*1.37 


2.62 
*2.67 


5.14 (-3) 
*5.20 


5.79 (-3) 
*5.82 


1.14 
*1. 14 


2.53 
*2.54 


5.0 


6.60 (-5) 
*6.59 


1.28 (-4) 
*1.27 


2.47 (-4) 
*1.94 


2.60 (-4) 
*2.59 


5.04 (-4) 
*5.02 


1.04 (-3) 
•1.03 


10.0 


*4.37 (-7) 


*8.26 (-7) 


♦1.64 (-6) 


*1.50 (-6) 


*2.83 C-6) 


*5.15 (-6) 


T 


6 = 1.0 
a=0.1 


6 = 1.0 

a=0.2 


6 = 1.0 
a =0.5 


6 = 1.0 
a = 1.0 


6 = 2.0 
a=0.1 


6 = 2.0 
a=0.2 




0.01 

0.02 

0.05 

0.1 


7.83 (-2) 

7.75 

7.66 

7.40 

7.00 


1.55 (-1) 

1.54 

1.52 

1.47 

1.39 


3.69 (-1) 

3.65 

3.60 

3.48 

3.29 


6.40 (-1) 

6.32 

6.24 

6.01 

5.65 


1.10 (-1) 
1.09 
1.08 
1.04 
9.70 (-2) 


2.19 (-1) 

2.17 

2.14 

2.06 

1.93 


0.2 


6.26 
*6.28 


1.24 
*1.24 


2.93 
*2.93 


4.99 
*5.0O 


8.53 
*8.53 


1.69 
*1.69 


0.5 


4.47 
*4.45 


8.85 (-2) 
*8.93 


2.07 
*2.07 


3.45 
*3.46 


5.82 
*5.79 


1.15 
*1. 15 


1.0 


2.56 
*2.57 


5.06 
*5.06 


1.17 
*1.17 


1.87 
*1.87 


3.14 
*3.14 


6.20 (-2) 
*6.21 


2.0 


8.49 (-3) 
*8.53 


1.67 
•1. 67 


3.74 (-2) 
*3.74 


5.61 (-2) 
*5.62 


9.57 (-3) 
*9.60 


1.88 
*1.88 


5.0 
10.0 


3.85 (-4) 
*3.29 

*1.67 (-6) 


7.55 (-4) 
*6.35 

*3. 16 (-6) 


1.77 (-3) 
*1.31 

*5. 77 (-6) 


7.99 (-3) 
*1.69 

*6.49 (-6) 


*3.37 (-4) 
•1.67 (-6) 


*6.54 (-4) 
*3. 16 (-6) 
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Table 4. Unscattered flux (4x/o-)I°(r,a,6) of radiation, at a corner position, from a rectangular -plane isotropic 
source calculated using tables 1 and 2 {numbers preceded by asterisks) and using table 3 (no asterisks) — Con. 





6 = 2.0 


6 = 2.0 


6 = 2.0 


6 = 5.0 


6 = 5.0 


6=5.0 


T 


a=0.5 


a— 1.0 


a = 2.0 


a =0.1 


a=0.2 


a=0.5 





5.25 (-1) 


9.31 (-1) 


1.41 (0) 


1.37 (-1) 


2.73 (-1) 


6.57 (-1) 


0.01 


5.18 


9.18 


1.38 


1.35 


2.68 


6.46 


0.02 


5.11 


9.05 


1.36 


1.33 


2.64 


6.35 


0.05 


4.91 


8.67 


1.30 


1.26 


2.51 


6.03 


0.1 


4.60 


8.07 


1.20 


1.16 


2.31 


5.55 


0.2 


4.03 


7.01 


1.02 


9.93 (-2) 


1.97 


4.72 




*4.02 


*6. 99 


*1.02 


*9.97 


*1.98 


*4.73 


0.5 


2.72 


4.62 


6.38 (-1) 


6.38 


1.27 


3. 00 




*2.72 


*4. 62 


*6.39 


*6.35 


*1.26 


*2.99 


1.0 


1.44 


2.35 


3. 02 


3.27 


6.46 (-2) 


1.50 




*1.44 


*2.35 


*3.02 


*3.27 


*6.45 


*1.50 


2.0 


4.24 (-2) 


6.43 (-2) 


7.46 (-2) 










*4.24 


*6.43 


*7.45 


*9.69 (-3) 


*1.90 


*4.28 (-2) 


5.0 
















*1.32 (-3) 


*1.75 (-3) 


*1.80 (-3) 


*3.37 (-4) 


*6. 54 (-4) 


*1.35 (-3) 


10.0 
















*5.79 (-6) 


*6.51 (-6) 


*6.52 (-6) 


*1.67 (-6) 


*3. 16 (-6) 


*5.79 (-6) 




6 = 5.0 


6=5.0 


6 = 5.0 


6 = 10.0 


6 = 10.0 


6 = 10.0 


r 


a = 1.0 


a = 2.0 


a=5.0 


a=0.1 


a=0.2 


a=0.5 





1. 19 (0) 


1. 88 (0) 


2.73 (0) 


1.47 (-1) 


2.92 (-1) 


7.06 (-1) 


0.01 


1.17 


1.84 


2.65 


1.44 


2.86 


6.92 


0.02 


1.15 


1.80 


2.58 


1.41 


2.81 


6.77 


0.05 


1.08 


1.69 


2.39 


1.33 


2.65 


6.38 


0.1 


9.92 (-1) 


1.53 


2.10 


1.21 


2.41 


5.79 


0.2 


8.35 


1.20 


1.64 


1.02 


2.02 


4.84 




*8.37 


*1.26 


*1. 65 


*1.02 


*2.02 


*4.84 


0.5 


5.14 


7.27 (-1) 


8.53 (-1) 


6.42 (-2) 


1.27 


3.01 




*5. 14 


*7.27 


*8.49 


*6.39 


*1. 27 


*3.01 


1.0 


2.47 


3.20 


3.44 










*2.47 


*3. 20 


*3.44 


*3. 27 


*6.47 (-2) 


*1.50 


2.0 
















*6.51 (-2) 


*7.57 (-2) 


*7.69 (-2) 


•9.69 (-3) 


•1.90 


*4.28 (-2) 


5.0 
















*1. 75 (-3 


*1.80 (-3) 


*1.80 (-3) 


*3.37 (-4) 


*6. 54 (-4) 


*1.35 (-3) 


10.0 
















*6. 51 (-6) 


*0.53 (-6) 


*6.54 (-6) 


•1.67 (-6) 


*3. 16 (-6) 


*5.79 (-6) 




6 = 10.0 


■ 6 = 10.0 


6 = 10.0 


6 = 10.0 


6 = 20.0 


6 = 20.0 


T 


a = 1.0 


a = 2.0 


a =5.0 


a = 10.0 


a=0.1 


a=0.2 





1.28 (0) 


2. 07 (0) 


3.15 (0) 


3. 80 (0) 


1.52 (-1) 


3.02 (-1) 


0.01 


1.26 


2. 02 


3.04 


3.64 


1.48 


2.95 


0.02 


1.23 


1.97 


2.94 


3.49 


1.45 


2.88 


0.05 


1.15 


1.83 


2.67 


3.10 


1.36 


2.70 


0.1 


1.04 


1.03 


2.30 


2.58 


1.22 


2.43 


0.2 


8.60 (-1) 


1.31 


1.74 


1.80 


1.02 


2.03 




*8.58 


*1.30 


*1.73 


*1.83 


•1.06 


*2.03 


0.5 


5.18 


7.34 (-1) 


8.65 (-1) 


8.78 (-1) 






1.0 


*5. 18 


*7.33 


*8.62 


*8.77 


*0.39 (-2) 


*1. 27 


*2.47 


*3.21 


*3.44 


*3. 45 


*3. 27 


*6.47 (-2) 


2.0 
















*6.51 (-2) 


*7.56 (-2) 


*7.69 (-2) 


*7.68 (-2) 


*9.69 (-3) 


•1. 90 


5.0 
















*1.75 (-3) 


*1.80 (-3) 


*1.80 (-3) 


*1.80 (-3) 


*3.37 (-4) 


*6.54 (-4) 


10. 
















*6.51 (-6) 


*6.53 (-6) 


*6.53 (-0) 


*6.53 (-6) 


*1. 67 (-6) 


*3. 16 (-6) 




6 = 20.0 


6 = 20.0 


b = 20.0 


6 = 20.0 


6 = 20.0 


6 = 20.0 


T 


a=0.5 


a = 1.0 


a=2.0 


a=5.0 


a = 10.0 


a = 20.0 





7.31 (-1) 


1. 33 (0) 


2. 17 (0) 


3. 38 (0) 


4. 22 (0) 


4. 88 (0) 


0.01 


7.13 


1.30 


2.10 


3.25 


4.01 


4.56 


0.02 


6.96 


1.27 


2.04 


3.12 


3.81 


4.28 


0.05 


6.50 


1.18 


1.88 


2.79 


3.30 


3.59 


0.1 


5.86 


1.05 


1.65 


2.36 


2.68 


2.81 


0.2 


4.86 


8.63 (-1) 


1.31 


1.75 


1.89 


1.92 


0.5 
1.0 


*4.85 


*8.62 


•1.31 


*1.74 


*1.86 


*1.89 


*3.01 


*5.18 


*7.34 


*8.63 (-1) 


*8.79 (-1) 


*8.82 (-1) 


•1. 50 


*2.47 


*3.21 


*3.44 


*3.45 


*3. 45 


2.0 
















*4.28 (-2) 


*6.51 (-2) 


*7.56 (-2) 


*7.69 (-2) 


*7.68 (-2) 


*7.68 (-2) 


5.0 
















*1.35 (-3) 


*1.75 (-3) 


*1.80 (-3) 


*1.80 (-3) 


*1.80 (-3) 


*1.80 (-3) 


10.0 
















*5.79 (-6) 


*6.51 (-6) 


*6.53 (-6) 


*6.53 (-6) 


*6.53 (-6) 


*6.53 (-6) 
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6. Examples 

Comparison of the results in the above sections 4 and 5 is made in table 5 with similar 
results of Berger and Lamkin [1] for a water barrier of t = 4 and a = b = 2, using the compound- 
ing rule in section 3 for noncorner positions. Agreement is within 5 percent. 

There exists little experimental information on flux fields from rectangular gamma-ray 
sources. However, an additional partial check on qo(a,b) was supplied by Auxier et al., [3] 
as shown by the circles, triangles, and squares in figure 7. The distributed source was rep- 
resented by an array of matched (±5%) Co 60 sources spaced 2 ft apart over a 24- by 32-ft rec- 
tangle. The detectors were pocket ionization chambers with an average response of 14.4 milli- 
roentgens per hour at 1 ft from a 1-millicurie point Co 60 source. This response calibration 
was made in a geometry such that the combined effects of backscattering from the ground, 
the buildup factor in air, and the exponential attenuation by air were roughly equivalent to 



Table 5. Dose rate, in milliroentgens per hour per milli- 
curies per square foot of source strength, at a number of 
positions a=x/h, £ = y/h at distance h from a rectangular 
surface of relative dimensions w/h= 2, l/h= 2 separated from 
a plane source of 1-Mev photons by a layer of water of 
thickness - = 4 

The numbers in parenthesis* are from table 1 of Berger and Lamkin [1]. 



\« 

\ 





0.25 


0.50 


0.75 


1.0 





0.211 
(.22) 

.275 

(.28) 

.323 

(.33) 

.350 
(.35) 

.359 

(.37) 










.25 


0.361 

(.37) 

.425 

(.44) 

.462 
(.48) 

.473 

(.49) 








.50 


0.503 
(.51) 

.546 
(.55) 

.559 

(.57) 






.75 


0.593 
(.60) 

.608 
(.62) 




1.0 


0.624 
(.63) 



*Arithmetic error in the Berger-Lamkin result at a— 0.25, 0=0.25 is here 
corrected. 




12 16 

FEET FROM CENTER, ALONG DIAGONAL 

Figure 7. Detector response in milliroentgens per hour at h = 3, 
5 and 11 ft above the diagonal of a rectangle of dimensions 
w = 24 ft, 1 = 32 ft covered with one millicurie per foot [2] of 
Co 60 . The circles, triangles, and squares are experimental 
values of Auxier et al. [3] and the curves are results from (28) 
compounded according to (11). 



these combined effects in the actual measurements. Thus, since the experimental data in 
figure 7 are related to a source strength of 1 millicurie per square foot, the dimensional factor 
<r/47r in (10) may be identified as 14.4 milliroentgens per square foot, and we may multiply 
qo(a,b; a,fi) by this number for an absolute comparison. The curves in figure 7 represent 
14.4-g(a,6; a,f$), and agree with the experimental data to within ±10 percent despite neglect 
of attenuation and scattering effects. 



The authors thank M. J. Berger and L. V. Spencer for their interest, encouragement, and 
helpful criticism regarding this effort. 

Appendix A. Evaluation of the p^ (a, b) Coefficients 

The Legendre polynomials, P z (cos 0), may be written in the form [13] 



fi(co S .) = g 2 ( !( ^; ( ^ ) !( co S .)-^ 



(34) 



where m = l/2 if I is even, and m=(l— 1)/2 if I is odd. Substituting (34) into (15) and using 
the notation 
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("> 



n\ 



(35) 

(n—r)\r\ 

for the r'th coefficient of an //'th degree binomial expansion, we have for the even p t 's 

„ 4,(-^pi-j\/4i-2j\r r» da dp 

lhl fy 2" V j An-j )Ja-oJi>-o(a t +r+iy->+*'*' 

Integration over a yields [14] 

' w (-!)«+' /2/-j\/4/-2j\/7-A „ +1 f» <& 



ft, -^&2*'(2t+l)i, j A2Z-iA i r Jo 



gi_j.l)i_i+i( a t + p» +1 )i+i 



(36) 



(37) 



For 1=0, integration over /3, using eq (234.19a,, reference [14], yields the result in (16). 
For the higher even terms, use of eq (234.18a) leads to (18) and higher terms which may be 
expressed in the general form 



Vu= 



r-iv-*-»W 2 * +2 *V l+k \ 

_L ab !=\ ( l) V l+k )\l-k-\) 

/(2/+l)2 2 'Va 2 +6 2 +lfe( 



■ (i) 

X(2l + 2k+l)j:Mf 



G>+ 



I) 



1 



U 2 2j (a 2 +b 2 +\y\(a 2 +iy- :i+i+ (b 2 +l) k -J +l )' 
The odd p t 's as obtained from (15) and (34) may be expressed as 

„ ^ (-iy (2i+i-j\/u+2-2j\ r r da dp 

Pu+x Z- 2 2 '+ 1 V j J\2l+l-jJJ a .oJ^o(a 2 +/3 2 +l)'^ +2 " 
Integration over a yields (see ref. 14 eqs 15, 15a, b) 

1 p tan" 1 (a/-y/p 2 +l)dp 



(38) 



(39) 



X 



l-j 

i-0 n! 



32 +1 )(-W/2 
1 



22i (2'-;^+l) (,_,-_<+!)£_, (^+IFT 



'/£ 



(a 2 +^ 2 +l)'-^ i+1 



(40) 



The integral containing the arctangent was amenable to integration by parts, and the purely 
algebraic integrand was reduced to terms of standard form by successive application of the 
identity 

(p 2 +l)(a 2 +p 2 +l) a 2 V/3 2 +l a 2 +,S 2 +l/ 
leading to the general form for the odd terms, including (17) and (19), 
1 



(41) 



P21 + V 



X 



_ ^ v-iy c fl\/l+k+l\ /2l+2k+l\ 
-2*'+\2l+l) fri, 2 2 * \kj\ k J\ l+k ) 



{1 b a i a 

{b 2 +\y v^+i tan_1 ^¥+i + W+W* V^+T tan ' 



Va 2 +1 



k 2 2j ~ 1 1/1 1 \1 

+ab ^j^)' (a 2 +b 2 +iy\la?+-Ly-J+ i+ (b*+iy-^ i ) J (42) 



Appendix B. Evaluation of the q n (a, b) Coefficients 

The integral (27) may be written in cartesian coordinates as 

I°(r,a,b) = ~ £da£dfi exp (-rVa 2 +£ 2 +l)/(c* 2 +/3 2 +l). (43) 

The substitution 

a=s/l3 2 +l tan 7] (44) 

yields 



tan - * - 



(45) 



/'(r,a,6)=f -^= exp(-rV/3 2 +lsecr;)^. 

4?r Jo V/3 +1 Jo 

Expanding the exponential as a power series of its argument we have 

The first integration is then simply over powers of sec rj, which can be found in integral tables 
as*standard forms [15], giving 

I°(r, a, b)-£ J o 6 { A W Un_1 W+i- T Sinh_1 W+l 

+2l a -^r ~2 "+ 2- Smh vTOJ+4! L3 +a( ^ +1) J-- • • •/*■ (47) 

Integration term by term over fi then yields a power series in r of the desired form (10), in 
which the coefficients g w (a,&) are 

q ° (a > h)= ^wn tan_1 wti **> (48) 

2l (a, 6)==- J q sinh-^== tZft (49) 

g 2 M)=|jJV etc, (50) 

as evaluated using (47) . 

No success was had in evaluating qo(a,b) in closed form. However, use of the arctangent 
expansion [16] 

, _ x - 2 2 "(m!) 2 / x 2 Y +1 , K1 . 

x tan x a;=S / Q i i m ( TT~2 ) ( 51 ) 

^r (2m + l)! \l+ry 

yields the series of integrable terms 

•*•«-! JTsKU$+i)" H » (52) 

Integration of these terms (see ref 15, p. 25) gives the result 



2o (a,&)=tan ^= • _ g ^ (^ 

+* -p 1 / « 2 Y +1 « 2*'(j!) 2 / a 2 +l V+ 1 

^a ito 1 2i+l VW ^ (2j+l)! W+b*+l) {06) 
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in which the double summation vanishes for i=0. Since convergence of the summations in 
(53) is poor when both a and b are large as compared with unity, it is advantageous to rewrite 
the first member of the right-hand side, using the substitutions 



tan VTO = 2- taQ 6" (M) 



and 

Va 2 +1 



2 ^TT ( TTT ) =tanh~ 1 ^sinh" 1 a. 

This gives 



(55) 



v« + l 



, ,. t . , , , _,Va 2 +l V« 2 +l^ 1 / a 2 V 

+ a&2i+lU 2 +l/ p{ (2j+l)\\a 2 +b*+lJ W 

which, after factoring the quantity 

&f, 1 /a 2 V +1 

from the last two members, may be identified with the result in eq (28). Now, since the inner 
summation over j in (28) consists of the leading terms of the expansion of (^a 2 +l/b) 
Xtan -1 (-yJa 2J rl/b) according to (51), the quantity within the braces in (28) is the remainder after 
i—\ terms of this expansion, causing the outer summation over i to converge very quickly for 
all values of a and b. 

The remaining q n (a,b) coefficients for n even and odd, respectively, may be written in 
the general forms 

1 f n ~ 2 /n — 1\ n- 1 /?) — 1\ a 2n-2i-l-k2i+\ -\ 

&(«, 6)=(^{(-l)- g (-D< ^ , >+2) W«, 6)+g (» . ') (2n _ 2i _ 1)(2i+1) } 

(57) 
and 

a. + i(a,ft)=-^ I j I {(2»-l)!^[-ft.-,(a,6)] 

+ W + i) [^ +1) * ■Wi-^-H^+i)- sioh- -^ 

+a6V« 2 +6 2 +l £ ^npjjj ( a 2 +6 2 +l)*{(a 2 +l)"- i - 1 +(6 2 +l)"- i - 1 }]| (58) 

from which g w (a,6)'s for 1 < n < 4 are given in (29 to 32). 
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